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LATTICE SUPERSYMMETRY WITH DOMAIN WALL FERMIONS 
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Supersymmetry, like Poincare symmetry, is softly broken at finite lattice spacing provided the 
gaugino mass term is strongly suppressed. Domain wall fermions provide the mechanism for 
suppressing this term by approximately imposing chiral symmetry on the lattice. We present 
the first numerical simulations of A/" = I supersymmetric SU(2) Yang-Mills on the lattice in 
d = 4 dimensions using domain wall fermions. 
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1 Introduction 

Supersymmetric (SUSY) field theories may play an important role in describing the physics 
beyond the Standard Model. Non-perturbative numerical studies of these theories could provide 
confirmation of existing analytical calculations and new insights on aspects of the theories not 
currently accessible to analytic methods, similar in spirit to the current situation in many other 
field theories, most notably QCD. Several, but not all, SUSY theories can be formulated and 
studied numerically on the lattice Bi'i. Foremost, as the lattice spacetime is discreet, only a 
discrete subgroup of the Poincare symmetry is unbroken, so any formulation will break SUSY. 
However, the allowed operators that break Poincare symmetry are irrelevant in the continuum 
limit, so one can calculate at several lattice spacings a and take the a — > limit. 

If a SUSY model has non-trivial scalar fields, scalar mass terms will break SUSY unless 
forbidden by some symmetry. Since these operators are relevant fine tuning will be needed to 
cancel their contributions in the continuum limit. We avoid this problem as the four-dimensional 
M = 1 Super Yang-Mills (SYM) theory does not involve scalars. The lattice fermion doubling 
problem will also break SUSY due to a mismatch in the number of bosonic and fermionic degrees 
of freedom. Removing doublers breaks chiral symmetry, allowing relevant gluino mass terms. 
In the traditional approach 0, fine tuning is used to cancel these mass termsu^ Pioneering work 
using these methods has already produced very interesting numerical results Q'c. 



Many properties of the = 1 SYM have already been computed using analytic techniques 
that take particular advantage of the supersymmetry in the model. On the contrary, numer- 
ical simulations of SYM are at least as computationally difficult as dynamical QCD. So, one 
must be careful to choose an interesting problem since the trade-off is to not study some aspect 
of dynamical QCD. For us, one such interesting problem is to determine if the vacuum sup- 
ports a non-zero gluino condensate (xx)) as widely believed, and whether the formation of the 
condensate is due to spontaneous or anomalous symmetry breaking. 

The Dirac operator in the adjoint representation of SU(A^) has an index 2Nu^ where v is the 
winding of the gauge field. Classical instantons have integer winding and cause condensations 
of operators with 2N gluinos which anomalously breaks the U(l) i?-symmetry to 7^2N- For 
(xx) to condense, the remaining 7j2n symmetry must further break either spontaneously or 
anomalously to Z2. If the breaking is anomalous, then the responsible gauge configurations 
must have fractional winding!. It has already been established that such gauge configurations 
do exist i. It is our goal to distinguish between these two scenarios. 

This work summarizes results recently presented by the author in collaboration with John B. 
Kogut and Pavlos M. VranasH. All numerical simulations were run on QCDSP supercomputers 
at Columbia Univ. and Ohio State. For reviews on DWF please see the LATTICE '00 review 
talk of VranasEi and references therein. The possible use of DWF in SUSY theories has been 
discussed in earlier works i'i and the methods used here are along these lines. For lists of 
references not included here for lack of space, please see the cited articles H. 

2 The DWF formulation 

In the lattice DWF formulation of a vector-like theory in d = 4 dimensions, the fermionic fields 
are defined on a d -|- 1 dimensional lattice using a local action. It is often convenient to treat 
the extra dimension as a new internal flavor space. Thus, the gauge fields are introduced in the 
standard way in the d dimensional spacetime and are coupled to the extra fermion degrees of 
freedom like extra flavors. The detailed form of the action can be found in our recent workQ. 

The key ingredient is the free boundary conditions imposed on the extra dimension. As a 
result, two chiral exponentially bound surface states appear on the boundaries (domain walls) 
with the plus chirality localized on one wall and the minus chirality on the other. The two 
chiralities mix only by an amount that is exponentially small in the size of the extra dimension, 
called Lg , and together form a Dirac spinor that propagates in the d dimensional spacetime with 
an exponentially small mass. Therefore, the chiral symmetry breaking artificially induced by 
the Wilson term can be controlled by the new parameter Ls- It is often convenient to introduce 
an additional bare gluino mass m/ to help control extrapolations to the chiral limit. In the 
Lg 00, m J ^0 limit chiral symmetry is exact at any lattice spacing without fine tuning. 

The computing requirement is linear in Lg, in contrast to traditional lattice fermion regula- 
tors where the chiral limit is approached only as the continuum limit is taken, a process that is 
achieved at a large computing cost. Specifically, because of algorithmic reasons, the computing 
cost to reduce the lattice spacing by a factor of two grows by a factor of 2^"^*^ in four dimensions. 
Therefore, the unique properties of DWF provide a way to bring under control the systematic 
chiral symmetry breaking effects using today's supercomputers. 

The application of DWF to A/" = 1 SU(iV) SYM is quite similar to iV/ = 1 SU(iV) QCD. 
The differences are merely that the fermions are in the adjoint representation of SU(A^) and 
that the Dirac fermion fields must satisfy the Majorana constraint, thereby reducing by half the 
number of fermion degrees of freedom. After integrating out the Grassmann fields subject to 
the Majorana constraint, the Dirac determinant is replaced by a Pfaffian: essentially the square 
root of the determinant, provided it is positive. For DWF, the Pfaffian cannot change sign so 
working with either Pfaffians or square roots of determinants is equivalent. 



3 Numerical results 



All numerical simulations were performed on 8^ and 4^ spacetime volumes using the inexact hy- 
brid molecular dynamics (HMD) R algorithm. The algorithm numerically integrates the classical 
equations of motion as part of generating a statistical ensemble with weights proportional to the 
fourth root of a two adjoint flavor Dirac determinant. For DWF, this weight is proportional to 
a single adjoint Majorana flavor. The step sizes were chosen such that systematic uncertainties 
due to numerical integration errors are negligible compared to statistical uncertainties. 



8", |i=2.3, M„=1.9 




Figure 1: Wilson line (\W\) on an 8'' lat- Figure 2: Gluino condensate (xx> vs. Ls 

tice. Diamonds are quenched and the cross is with /3 = 2.3 and m/ = 0.02, 0.04, 0.06, 0.08. 

dynamical with Ls = 24 and m/ = 0. Curves are an exponential fit. 



The 8^ volume simulations were done with /? = i/g^ = 2.3. The value of /3 was chosen as 
large as possible without entering the thermal transition region. From quenched simulations 
in an 8^ volume, the magnitude of the fundamental Wilson line the quenched order 

parameter, is plotted vs. (3 in figure |^. The value of from a simulation of the dynamical 

theory at /? = 2.3 is also shown (cross), indicating that the dynamical theory is in a phase that 
"confines" fundamental sources. Therefore, the box size is large enough to avoid thermal effects 
that break SUSY. Using similar arguments, the 4^ volume simulations were done with /? = 2.1, 
near the limit of the weak coupling regime. Scaling arguments appropriate for weak coupling 
suggest that the lattice spacing at /? = 2.1 is twice as large as at /? = 2.3. 

To extrapolate the measured values of {xx) to the chiral limit, Lg ^ oo and ^ 0, 
simulations were performed in 8^ volumes at fixed (3 = 2.3 while the size of the extra dimension 
Ls was varied between 12 and 24 and the bare mass mf was varied between 0.02 and 0.08. 
The measured values appear as the points in figure |^. According to Leutwyler and SmilgaEI, 
if the formation of the gluino condensate is due to spontaneous symmetry breaking, the lattice 
volume limits how small the dynamical qluino mass can be set without losing the condensate: 
12meff (xx) V ^ 1 (the 12 is just normalization). As m^s ^ mf, this limit is satisfied for all 8^ 
simulations with ruf > 0.02. 

To estimate the gluino condensate in the chiral limit, we first extrapolate at fixed mj to the 
Lg oo limit using the fit function {xx) = co + ci exp(— 02^^). The best fit functions are plotted 
as the curves in figure |2[ The values of the extrapolated gluino condensate (with propagated 
errors) appear as points in figure ^. These extrapolated values are then further extraploted to 
the mj — > limit using a linear function (xx) (-^s — > 00) = 69 + 6imj. The best fit function 
appears as the line in figure ^. It is also reassuring to note that reversing the order of limits, 
i. e. first ^ at fixed Lg then Lg — > 00, yields a statistically consistent answer. 

Another approach to estimating the gluino condensate in the chiral limit is to actually 
perform dynamical simulations with m/ = 0. Since finite Lg will induce an exponentially small 
breaking of chiral symmetry, the effective gluino mass will not be zero. However, the gluino 
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Figure 3: Extrapolated (xx) to is ^ oo 
limit vs. nif and linear fit to m/ ^ limit. 



Figure 4: Dynamical (xx) at nif — vs. Ls 
on an 8* lattice at P — 2.3 and a 4* lattice 
with (3 — 2.1. Curves are exponential fits. 



mass should be too small to support spontaneous symmetry breaking. Additional simulations 
were run for Lg = 12, 16, 20, 24. The data are shown in figure The curve is the best fit to the 
exponential fitting function with the extrapolated value of the condensate as shown. 

Surprisingly, both methods for estimating the gluino condensate produce consistent results 
within the statistical errors. Note that this is inconsistent with the notion of spontaneous 
symmetry breaking. Operationally, this result reinforces our claim that systematic uncertainties 
are still relatively small despite limited statistical precision. Further, it gives us some confidence 
that our fit functions are valid over the region of interest. 

To further check for spontaneous symmetry breaking of the Z4 symmetry, we measured (xx) 
on smaller 4^ lattices with mj = and even larger values for Lg- The data are shown in figure]^ 
with the best exponential fit and the extrapolated value for the condensate. This provides even 
stronger evidence that spontaneous symmetry breaking is not responsible for the formation of 
a gluino condensate, at least in finite volumes. On these lattices 12m fV {xx) < li so analytical 
considerations!! suggest the support of (xx) must come primarily from topological sectors with 
fractional winding of 1/ = ±1/2. 

The spectrum of the theory is of great interest but it was not possible to measure it on the 
small lattices considered here. Also, the gluino condensate was measured at only two different 
lattice spacings so extrapolation to the continuum limit to compare with analytical results is 
not possible. Future work could explore these very interesting topics. 
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